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Abstract 

This is a pedagogical account of the some of the recent results of Brydges 
and Imbrie, described from the point of view of Grassmann integration. Some 
simple extensions are pointed out. 

Historical Introduction 

The subject of dimensional reduction, in the context considered in this lecture, has a 
long and rather peculiar history. The strangest aspect of all is that the simplest formulation 
of the results has been found only very recently. 

That one physical theory in d space dimensions should be related to another m D = d — 2 
dimensions was first suggested in the context of the field-theoretic formulation of the critical 
behaviour of the Ising model in a quenched random magnetic field [1,2]. It was noticed that 
the most singular Feynman diagrams all had the form of trees before averaging over the 
quenched randomness. When this was done, it turned out that the diagrams were the same 
as those for the critical Ising model without any random field, in two fewer dimensions. 
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Thus, no further work was needed - one should be able to take, for example, the known 
results for the ordinary Ising model in D dimensions, and find the critical behaviour of the 
random field Ising model m. d = D + 2 dimensions. In particular, since the lower critical 
dimension (at and below which there is no phase transition) of the ordinary Ising model 
is D — 1, this should imply that the 3-dimensional random field Ising model has no phase 
transition. 

Unfortunately this beautiful idea is wrong. In particular it contradicts a very simple 
argument of Imry and Ma [1] that the lower critical dimension should be two. However, it 
took several years of theoretical and experimental confusion and contradictory papers before 
Imbrie [3] proved that Imry and Ma were right. 

During this time Parisi and Sourlas [4] came up with a startling explanation of how 
and why dimensional reduction works: supersymmetry. At that time supersymmetry was 
thought to be the preserve of particle physics and string theory, so their argument was not 
widely understood. In their paper was also the germ of the reason why simple dimensional 
reduction might not work for the random field Ising model: the tree diagrams were seen to 
be the perturbative solution to the classical field equations in the presence of the random 
external field. If the solution of these equations was unique and given by the sum of the 
diagrams, then all should have been well. But it was easy to see that, at low enough 
temperatures, the solution was not unique, and that the one with the lowest energy was not 
the perturbative one. Work is still continuing on how, possibly, to rectify this situation, but 
this seems a very difficult problem. 

Meanwhile, a few years later, Parisi and Sourlas [5] came up with another ingenious 
application of their ideas. The problem known under various versions as branched polymers, 
lattice trees, or lattice animals (to be decribed later) had been attracting some attention. 
Lubensky and Isaacson [6] had formulated a rather complicated field theory for this problem, 
had come to the conclusion that the upper critical dimension (above which mean field theory 
is valid) is li = 8, and had noticed that the first terms in the e-expansion of the critical 
exponents below this were the same as those of the so-called Yang-Lee edge singularity 
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in D — d — 2 dimensions. This is the problem of an Ising model in a purely imaginary 
magnetic field, and is described by a simple scalar field theory with a cubic interaction 
and a purely imaginary coupling. Parisi and Sourlas reformulated the model of Lubensky 
and Isaacson in a simpler way, and then used their supersymmetry to explain why the 
dimensional reduction happened. This time, everything was all right. The numerically 
measured exponents of lattice animals and lattice trees in c? = 3 dimensions are given by 
the exactly known Yang-Lee exponents in D — 1. A non-perturbative proof of dimensional 
reduction in a supersymmetric field theory was given [7]. It was understood that the main 
reason for the failure of dimensional reduction in the random field Ising model does not 
apply here, essentially because the tree diagrams are the branched polymer configurations 
[8]. But none of these arguments were completely watertight, largely because it seemed 
necessary to ignore so-called irrelevant terms which spoiled the supersymmetry, but were 
believed not to affect the critical behaviour. 

In fact, the Yang-Lee critical theory is nowadays strongly beheved [9] to be in the same 
universality class as an even simpler problem: a classical gas with short-range repulsive 
interactions in the grand canonical ensemble. As a function of the fugacity 2; > 0, the grand 
partition function is sum of positive terms and has no singularity nor any zeroes. But in the 
complex z plane the closest singularity to the origin lies on the negative axis ai z = —Zc, 
and the critical behaviour near this is believed to be independent of the particular details 
of the interaction (so it is known as the 'universal repulsive gas singularity'), and it is also 
believed to be the same as that of the Yang-Lee problem at the critical imaginary magnetic 
field. Thus the branched polymer problem in d dimensions was believed to be related to the 
universal repulsive gas singularity in D — d — 2 dimensions. But once again these arguments 
relied on the neglect of irrelevant terms in the appropriate field theories. 

It took Brydges and Imbrie [10] to realise that most of this theoretical undergrowth could 
be cut away - that field theory is not needed at all! There is a simple and physical model 
for branched polymers in d dimensions for which dimensional reduction to a repulsive gas 
is rigorous and exact, not just in the critical region (which in this context means very large 
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polymers) but for all values of the fugacity. The arguments do use supersymmetry, in a very 
beautiful way, and give compelling evidence of how important it is in theoretical physics to 
learn how to walk (do problems with finitely many particles) before one tries to run (field 
theory and beyond.) 

The aim of these lectures is to present these ideas in as simple a manner as possible, 
using the ideas of Grassmann integration nowadays familiar to many theoretical physicists. 
I shall, however, quickly recall all the ingredients necessary for the argument. 

Classical gas and cluster expansion. 

Consider a classical gas of point particles, in some large D-dimensional domain fl, inter- 
acting through a rotationally symmetric two-body potential V. The grand partition function 
is 

oo 

E{z) = E Znz'' (1) 
Ar=o 

where 




where = {rj — r^)^ and we work in units where kT = 1. As the volume of the domain 
goes to infinity, we expect, if V is sufficiently short-ranged, that S(^) ~ e^*-^'"^', where p{z) 
is the pressure. 

For our purposes it is easier to consider the mean density n = z{d/dz)p{z), which has 
the expansion 

oo N I- N 

n{z) = {Eiz))-' Y: Jl^rui L n d^r,6^{n)e-^^<^''^'^^^ (3) 

N=l \ )■ ^ j=l 

The cluster expansion follows by writing e~^^^^^ — 1 — G{r^), so that e~^J<''^^^^'^^ — 
Ylj^k (1 ~ ^('^jfe))) expanding everything in powers of G. Note that if V is short-ranged, 
so is G. For the density, all the contributions which diverge as \fl\ — > oo cancel between 
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the numerator and denominator, and we are left with a sum of connected cluster diagrams. 
The coefficient of is a sum of connected graphs Q, each of which has N nodes labelled 
by (ri, . . . , rjv), of which ri is fixed to be at the origin. The value of the graph is given by 
integrating over all the other Vj the product of the —G{r'j^) over all connected edges (jk) of 
the graph. An example is shown in Fig. 1. 

Note that the integral corresponding to a particular graph Q can be performed by in- 
serting the integral representation 

/•CO 

G{r%) = / /(a,,)e-"^-^'=da,, (4) 

J 

on each edge: the integral over the Vj then has the form 
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where Ag is a matrix which depends on the {djk}-, and on but not on D. The only feature 
of this we shall need is that it factorises into D integrals over each component of the vectors 
Tj, and by the rules of gaussian integration it has the form 

(27r)(^-W(detAg)-^/2 (6) 

Grassmann integration 

Grassmann coordinates behave in many ways like ordinary coordinates, except that they 
anticommute. For our purposes, we are going to consider a space with just two Grassmann 
coordinates, called 9 and 0. It is useful to consider these as analogues of complex commuting 
coordinates z — x -\- iy and 'z — x — iy, so that the squared distance between two points 
{6j, 9j) and {9k, 9k) is the analogue of {zj — Zk){zj — Zk), that is \9jk\'^ = {9j — 9k){9j — 9k). 

We need the notion of integrating some function f{9, 9) over all space. Note that because 

2 

9 — 9'^ — 0, the Taylor expansion of any such function terminates: f — A + B9 + C9 + D99. 
One important property of integration over all space is that we may shift the integration 
variables by constants, so that 
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J d9d9f{e + K,e + K)^ J dedef{e, e) (7) 

This can only be satisfied if 

j dm 1 = j dm 6^ j dm e = (8) 

but we are allowed to choose the normalisation of 

jdmee^K (9) 

We also have to learn about gaussian Grassmann integration: the main rule is that if B 
is an N X N matrix 

J n dejddje-'^ EjJjB,,ek ^ (_1A)^ det B (10) 
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We shall not try to prove this in general [11]. Let us note the simplest case N — 1: 

J dddee-^^^^ = J dede{i - \Bee) = -\kb (ii) 

Supersymmetric classical gas 

We are now going to consider a classical gas which lives in superspace. Points in super- 
space are labelled by d commuting coordinates, the components of r, and two anticommuting 
coordinates {9, 9) as above. The squared distance in superspace is 

R]k = r% + \9,k? = (r, - nf + {9, - h) {9, - 9^) (12) 

and our classical gas with N particles has grand partition function Sss(-z) = S?/=o Zn,ssz^ , 
where 

Zn,ss = J^ JII d'r,d9,d9,e- ^^<^ ^^^^'^^ (13) 

Here 'SS' stands for superspace, or for 'supersymmetric'. The supersymmetry in this case 
is under super-rotations in the d + 2 dimensional space. Note that there are many different 
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kinds of supersymmetry. This one is rather different from that considered in most particle 
physics apphcations, and is also different from the global super symmetries used in quenched 
random systems in condensed matter physics. 

Actually, the partition function for this supersymmetric gas is not very interesting: since 
V depends only on the relative coordinates, the integration over the center-of-mass anticom- 
muting coordinates gives zero, by virtue of (8), except when A'" = 0. Thus Sss(;z) = 1. This 
is typical of supersymmetric theories. However, the mean density is non-trivial: 

N=l j=l 

where 5{9i)S{9i) just sets = = in whatever expression follows it. 

We are going to evaluate (14) in two different ways: first we shall show that it is identical 
to the same expression (3) for the ordinary classical gas in D = d — 2 dimensions, then we 
shall evaluate it in another way to show that it gives branched polymers in d dimensions. 



Dimensional reduction 

As before, we can make a cluster expansion of (14) by writing = 1 — G. In this 
case, only connected graphs survive because any disconnected pieces vanish on integrating 
over their center-of-mass anticommuting coordinates. Using the integral representation (4), 
a given graph Q now contributes (before doing the {ajk} integrations) 

/ l[d''rjdejdej6^''\ri)6{ei)6{ei)e-^^^^ (r,A,,r,+e,A,,e,) 
= (27r)(^-iW2 (det Ag)-'^/2 {-\KY-\detAg) 

and we see that this is identical to (6), with D — d — 2, as long as we take 

A = -l/7r . (15) 

To summarise, we have 

4f)(^) =n('^-2)(^) (16) 
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In a sense, the integrations over the anticommuting coordinates cancel those over two of the 
commuting ones. This again is a result typical of supersymmetry. 

Branched polymers 

We haven't yet tried to define the term 'branched polymer.' Physically, it is an object 
consisting of nodes connected by linear segments which are chains of monomers. Each node 
has a specific degree, the number of segments meeting at that node, which can either be 1 
or any integer > 3. By extending this defintion, we can think of the connections between 
each pair of monomers on the linear sections as nodes of degree 2. Each branched polymer 
has the topology of a tree, with a fixed number of monomers on each segment or branch. 
It is embedded in d-dimensional euclidean space, and there are assumed to be interactions 
between the monomers which (a) represent the fact that neighbouring monomers on the tree 
are confined within a typical radius a, and (b) there is a steric repulsion between all other 
pairs of monomers which inhibits them from approaching each other closer than a distance 
also 0(a). The branched polymer is assumed to be in thermal equilibrium with some kind 
of solvent, so that is explores its allowed phase space, each configuration occurring with a 
probability oc Q-^^^^sy/kT^ according to the laws of statistical mechanics. 

For enumeration purposes it is often simpler to allow the nodes to lie at the sites of a 
regular lattice, of spacing 0{a), and then to count, with equal weights, all allowed configu- 
rations of a given total number of monomers N. Depending on whether loops are allowed 
or not, this gives the problems of lattice animals or lattice trees. Numerical work sug- 
gests that, for large A^, these models exhibit a kind of critical behaviour, characterised by 
universal critical exponents. For example, the typical size R ~ N", and the number of dis- 
tinct animals (or trees) which are not related by lattice translations pn, goes like N'^~^jji^ , 
where ii is model-dependent, but the exponents v and a are universal. The connection with 
conventional critical phenomena is made more clear by considering the generating function 
g{^) — J2nPnz^ , where z — Zc — IJ>~^ is the critical point, near which g{z) ~ {zc — z)^'"". 
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We shall show that the supersymmetric model (13) gives rise to a continuum model of 
branched polymers which in fact is much more realistic than these lattice models. 



From the supersymmetric gas to branched polymers 

Go back to the formula (14) for the density and write, for each pair of particles 

e-^(«.-) ^ P{R%) = P{t% + = P{r%) + g(rj,)^,,^,, (17) 

where (5(r^) — P'ir"^). So we can write the integrand in (14) as a sum of 2^(^~^)/^ terms, in 
each of which for each pair [jk) we choose either the term P{r'jf,) or the term Q{r'j^9jk9jk. 
Each term can be represented by a subgraph of the complete graph connecting all N nodes, 
in which the edge {jk) is included if we choose the second term, and excluded if we choose 
the first. See Fig. 2. Now consider performing the integrations over the anticommuting 
coordinates. For the same reason as above, any subgraphs not connected to the origin will 
vanish. Moreover, it is fairly easy to see that any subgraphs containing cycles will also vanish, 
since this will always involve some anticommuting coordinate being raised to a power > 2. 
Thus the surviving subgraphs all have the form of rooted connected trees, rooted because the 
particle at the origin is privileged. It also easy to see that all possible such rooted trees with 
N nodes occur (we allow nodes of degree 2.) For a given tree T we may then relabel the 
coordinates, starting from the origin. The integration over the anticommuting coordinates 
on one particular edge of the subgraph is then of the form 

j dedeQ{r^)ee = {-l/7r)Q{r^) (18) 

We therefore have 

oo 

4f ^(^) = E zi-z/Trf-'NZ^ (19) 

N=l 

where 

1 r ^ 

^AT = E ICn / n d'r.S^'^ (ri) n Q(r|,) U P{r%) (20) 
r •' j=i j~fe j/fc 
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where j ^ k means that j and k are neighbouring nodes on T, and j k means that they 
are not, and each pair {jk) is counted just once. 

For suitable potentials V, we can interpret this as the partition function for a branched 
polymer with the topology of T: P{r'j^ — e"^^''^*^ is the Boltzmann weight representing the 
steric repulsion between different parts of the polymer, and Qir"^/.) then is the weight for 
two neighbouring nodes, which keeps them a distance 0(a) apart. If the original gas had 
only repulsive forces, then P' > so that Q may be interpreted as a weight. Fortunately 
it is quite easy to arrange a suitable function Q. An example is illustrated in Fig. 2. A 
special case is to choose V{r'^) = Vq > for r < a and zero for r > a. Then 

Q{r^) = (1 - e-^")5(r' - a") = {2a)-\l - e-^')6{r - a) (21) 

If we take the limit Vq — > oo corresponding to a hard core repulsive gas, then the correspond- 
ing branched polymer model may be thought of as consisting of balls of radius a/2, with 
neighbouring balls on the tree constrained to touch each other, and a hard core repulsion 
between the rest. 

Putting together all the pieces of the puzzle, we have shown that the generating function 
for rooted branched polymers in d dimensions is related to the mean density of a repulsive 
gas in D — d — 2 dimensions by 

oo 

z{d/dz)p{z) = n{z) = z{-z/t:)^-^NZn (22) 

N=l 

By integrating with respect to we finally end up with a relation between the pressure 
p{z) — lim|n|^oo If^l"^ lnS(2;) of the classical gas and the generating function for branched 
polymers: 

p{z) = -7rE(-z/7r) (23) 

where 

oo 

E(C) = E ZatC"^ (24) 

Af=l 

Eq. (23) is the main result of Brydges and Imbrie [10]. 
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Simple examples 



1. d = 2 

When d — 2 then D = 0, a particularly simple case, since the universe is all at one point 
[12]. In that case, wc do not have to worry about the thermodynamic limit either. The 
grand partition function is 

S(-) = E ^e-^^(^-)^(°) (25) 

since if there are N particles at the origin, there are \N{N — 1) interactions between them. 
The remarkable thing is that the branched polymer partition function in c? = 2 depends only 
on the value of the potential V{r) at the origin. This is a consequence of the supersymmetry. 
If we relax the condition that Q — P' , it is no longer true. Although the sum in (25) cannot 
be performed exphcitly, it can be shown [13] that all the zeroes of r.{z) he on the negative 
real axis and are distinct. The closest one to the origin determines the large N behaviour 
of 2^jv, and gives rise to the branched polymer singularity. 

In the hard core case V{Q) — > oo, there can only be zero or one particles, so we have 

p{z) = In S(z) = ln(l + z) (26) 

so that 

E(C) = -(l/7r)ln(l-7rC) (27) 

and^ 

= TT^-ViV (28) 

It is very hard to see how this simple result could be obtained by other means. While it is 
easy to check for = 1, 2, it already is non-trivial for N = while for iV > 4 the integrals 
involved in each contribution appear intractable. 



^Note that this differs by factors of 2 from equation (1.4) of Ref. [10], since we have Q = (l/2a)5(r— 
a), while these authors take Q = 8{r — I). 
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2. d = 3 



Interacting classical gases in D — 1 can generally be solved exactly only in the case of 
infinite hard core repulsion, i.e. V{0) oo. We can think of the particles as rods which 
we now take to be of unit length. If the total length is L and we use periodic boundary 
conditions then since configuration space has volume {L — N)^,we have — {L — N)^/N\ 
A classic calculation in statistical mechanics then shows that, in the thermodynamic hmit, 
p — n/ (1— TT-) (the denominator is just the excluded volume term in van der Waals' equation!) 
and so z{dp/dz) — n — p/{l+p). Integrating, we find an implicit equation for the pressure: 
z = pe^. Remarkably, this can be solved explicitly to find p{z) as a power series in z [Hint: 
write § p{z)dz/ z^^^ as a similar contour integral over p] 

Pi^) = E -A^^"" (29) 

N=l ■ 

Using the main result (23) this expansion once again gives some highly nontrivial identities 
for sums of certain integrals in three dimensions. 

If we are just interested in the branched polymer singularity note that z{p) has a 
quadratic turning point at z = —Zc = — e~^, so that p{z) has a square root branch point 
there. This is consistent with the square root singularity at the Yang-Lee edge in D — 1: 
it is an easy exercise in transfer matrices to solve the nearest neighbour Ising model in a 
purely imaginary field, and check this. 

Correlation functions 

Although we have so far considered only the mean density, it is easy to generalise di- 
mensional reduction to correlation functions: 

(n(ri)n(r2) . . . n(rp))g') = {n{n)n{r2) . . . n(rp))('^-''°) (30) 

which is valid as long as the points rj all lie in the {d — 2)-dimensional subspace, and at 
9 — 9 — 0. If we now expand the left hand side in terms of P and Q as before we find the 
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generating function for p different branched polymers, which exclude each other's volume, 
rooted at the points rj. Of course, for p 2 the correlation function enjoys full rotational 
invariance in superspace. Thus if we take ri = and ri = r and we let 

G{r^ + QQ) = (n(O)n(r, 6, (31) 

then G{r'^) is the same as the density-density correlation function for the D-dimensional 
ordinary gas. Writing now G(r^ + 99) — G{r'^) + 99G'{r'^), the first term gives as before 
the correlation function between two different but interacting branched polymers, while, 
after performing the Grassmann integrations, we see that the second term, proportional to 
G'(r^), gives the density-density correlation function on a single branched polymer. That 
these are related is a consequence of the supersymmetry. This is all spelled out in full detail 
in Ref. [14]. Since for D = 1 the correlation function can in some cases be computed exactly, 
so can that for branched polymers in o? = 3. This was first carried out using the Yang-Lee 
model by Miller [15]. 

Extensions 

I'll finish by briefly mentioning some simple extensions of this beautiful idea. 

Bells and whistles 

It is possible to decorate the nodes of the branched polymer with other degrees of freedom, 
e.g. to consider that they are of different species, and that the interactions between them 
can depend on this. Under dimensional reduction this becomes a classical gas with the 
same decoration. It would be nice if we could find a different universality class at negative 
fugacity. Field theory suggests that this might be the case (we choose a theory with a 0^'*+^ 
interaction with n > 2, rather than 0^) but investigations so far [13] have shown that, at 
least for physical values of the interactions, the universal repulsive singularity is robust, and 
nothing new happens. 
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Adsorption at a wall or a plane 



The astute listener should have reahsed that all that is needed for dimensional reduction 
to work is rotational invariance in two of the d commuting dimensions, so that integrals 
over these can be cancelled by the two anticommuting ones. Thus, for o? > 3, we are free 
to add an arbitrary potential U{x) depending on one of the commuting coordinates x. The 
generating function for branched polymers in d dimensions in the presence of this potential 
will be related to that of a D-dimensional classical gas in the same manner as (23). 

As an example take a hard core repulsive gas in d — 3, and e~^^^^ = 1 + X6{x), with 
A > 0, representing an attractive potential on the plane x — 0. A simple calculation then 
shows that the density at a point close to, but not on, the plane is modified to 



where no(^) is the unperturbed density. For small A, the nearest singularity to the origin 
is still at z — —Zc, the d — 3 branched polymer singularity. However, for large enough A it 
instead comes from the vanishing of the denominator, and is generally a simple pole. This is 
the branched polymer singularity in c? = 2. In this regime the branched polymer is confined 
to the vicinity of the plane x — 0. The critical value of A at which these singularities collide 
represents an adsorption transition. The case of an attractive wall, where U — +oo for 
X < 0, may be dealt with similarly. 



We do not have to stop at adding one pair of anticommuting coordinates - for example 
we can add another, in which case we get dimensional reduction by four. But the equivalent 
d-dimensional model then involves interaction weights proportional to P, P' and P", and 
the last typically has to change sign, so is unphysical. 




(32) 



Dimensional reduction by 4 
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Intriguing Connections 

If one takes the scaling limit of the expression (25) as z ^ —z^ and y(0) — > 0, one finds 
something proportional to an Airy function. The Airy integral arises in many contexts, 
but another was recently found by Richard, Guttmann and Jensen [16], who argued that 
it occurs in the scaling limit of self-avoiding closed loops in the plane, weighted by their 
length and the area enclosed within. The parameters conjugate to these are the fugacity, 
analogous to z, and the (negative) pressure inside the loop, which plays the same role as 
y(0) in the branched polymer generating function. Using the field theory arguments of 
Parisi and Sour las [5], I have argued elsewhere [17] why an Airy function is expected in this 
context. This also explains why the fractal dimension of self-avoiding loops is |. But it 
would be very interesting to provide a direct connection between branched polymers and 
self-avoiding loops - is perhaps the external perimeter (suitably defined) of a large branched 
polymer in the same universality class as that of ordinary self-avoiding loops? 
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FIGURES 




FIG. 1. A typical cluster diagram for the mean density. 




FIG. 2. A branched polymer configuration. Solid lines correspond to the nerest neighbour 
interaction Q, and dashed lines to repulsive interactions P. 




FIG. 3. Typical forms for P (dashed curve) and Q = P' (solid). 
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